Inter-layer Josephson coupling in stripe-ordered superconducting cuprates 
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Motivated by experiments on La2-ccBa x Cu04 which suggest that stripe order co-exists with two- 
dimensional pairing without inter-layer phase coherence over an extended range of temperatures, 
we determine the inter-layer Josephson coupling in the presence of stripe order. We employ a 
mean-field description of bond-centered stripes, with a zero-momentum superconducting condensate 
and alternating stripe directions pinned by the low-temperature tetragonal (LTT) lattice structure. 
We show that the Fermi-surface reconstruction arising from strong stripe order can suppress the 
Josephson coupling between adjacent layers by more than an order of magnitude. 



I. INTRODUCTION 



Stripe order is a fascinating phenomenon in cuprate 
superconductors i* 2 - Originally detected in neutron- 
scattering experiments on La 2 _ x _ 2y Nd 2y Sr :c Cu04, 3 this 
combination of uni-directional spin and charge order was 
found in other members of the "214" family of cuprates 
as well. 2 Remarkably, incommensurate low-energy spin 
fluctuations, often interpreted as precursors to stripe or- 
der, are seen not only in La 2 _ x Sr :c Cu04, but also in 
YBa 2 Cu 3 06 + 5 415 and Bi 2 Sr 2 CaCu 2 8+(5 i£ Together with 
STM measurements, which detected signatures of charge 
stripes (albeit with substantial disorder) on the surface of 
Bi 2 Sr 2 CaCu 2 8 +<5 and Ca 2 _ x Na :E Cu0 2 Cl 2 ^ these find- 
ings suggest that the tendency toward stripe order is com- 
mon to underdoped cuprates. 

What is less clear is the role of stripes for superconduc- 
tivity. A large body of experiments appears consistent 
with the concept of competing superconducting and mag- 
netic order parameters^ 2 - including e.g. the magnetic-field 
enhancement of spin-density wave (SDW) order. How- 
ever, a few observations also point to a co-operative in- 
terplay of SDW and pairing. In La 15 / 8 Ba 1 / 8 Cu04, the 
onset of SDW order upon decreasing temperature is ac- 
companied by a significant drop in the in-plane resis- 
tivity, while bulk Meissner effect sets in at much lower 
temperatures^ This intermediate-temperature regime 
of Lai5/ 8 Ba 1 / 8 Cu04 has been interpreted in terms of 
fluctuating 2d pairing, without inter-layer phase coher- 
ence. (A related phenomenon is the suppression of 
the Josephson plasma resonance seen in the optical- 
conductivity measurements of La 2 _ x Sr :c Cu04 upon ap- 
plication of a moderate c-axis field. 10 ) In order to explain 
the absence of an effective inter-layer coupling, the exis- 
tence of a stripe-modulated (i.e. finite-momentum) su- 
perconducting condensate, a so-called pair density wave 
(PDW), was postulated ] 11 ! 12 Indeed, in the absence of a 
uniform condensate, the lowest-order Josephson coupling 
between neighboring layers vanishes, if the condensate 
modulation direction alternates from layer to layer - the 
latter being the result of the in-plane lattice distortions 
inherent to the LTT structure. However, some properties 
of the PDW state are not easily compatible with exper- 
iments: both photoemission and STM have established 



a (i-wave like gap in the stripe-ordered state above T c , 
whereas the PDW state has a full Fermi surface. 

In this paper, we shall investigate whether a primar- 
ily uniform condensate in a stripe-ordered state could 
be compatible with the scenario of fluctuating 2d pair- 
ing in La 15 / 8 Ba 1 / 8 Cu04, i.e., the absence of inter-layer 
phase coherence. To this end, we calculate the inter-layer 
Josephson coupling for a mean-field model of supercon- 
ducting charge-density wave (CDW) state with realistic 
parameter values. The stripe order induces a reconstruc- 
tion of the Fermi surface, with rotation symmetry break- 
ing in each Cu0 2 plane. As the stripe direction alternates 
from layer to layer, the reconstructed Fermi surfaces of 
adjacent layers do not match in momentum space. This 
effect leads to a significant suppression of the Josephson 
coupling; an additional suppression arises from incom- 
mensurate antiferromagnetism of realistic amplitude. 



II. ORDER PARAMETERS AND MEAN-FIELD 
MODEL 



We start by enumerating the relevant order parameters 
for a superconducting stripe state. Charge and spin den- 
sity waves, with wave vectors Q c and Q s , are related to 
expectation values of particle-hole bilinears in the singlet 
and triplet channel, respectively. For instance, a CDW 
is characterized by non-zero F c (k) = Xlcr( c | + Q a c ko) 
where ct creates an electron with momentum k and 

kcr 

spin a. The superconducting condensate can have both 
a uniform component and a modulated (PDW) compo- 
nent with wavevector Q p , such that (c£_|_q ^ c %\) ^ 0- 
On symmetry grounds, a collinear SDW will induce a 
CDW with wavevector Q c = 2Q S , a PDW will induce a 
CDW with Q c = 2Q P , and in the presence of a uniform 
condensate a CDW will induce a PDW with Q p = Q c - 

In our modelling, we start from two Cu0 2 layers with 
homogeneous d-wave pairing and then add CDW/SDW 
modulations. Each layer i = 1, 2 is described by a quasi- 
particle model of electrons moving on a square lattice, 
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FIG. 1: Schematic real-space structure of valence-bond 
stripes at doping 1/8,— with circle sizes (line widths) denot- 
ing on-site hole densities (bond strengths), a) Paramagnetic 
state with dominant d-wave modulation, b) Additional spin 
modulation with anti-phase domain walls on hole-rich stripes. 



with the Hamiltonian 



(i) 



and the in-plane dispersion = —2t(cosk x + cosk y ) — 
At' cos k x cos fc y — 2t"(cos2k x + cos 2/^) with t' = — 1/4 
and t" = t/12. 

The symmetry-breaking orders are implemented non- 
selfconsistently at the mean- field level by Hdw and Hp . 
We restrict our attention to bond-centered stripes of 
period 4 (8) in the charge (spin) sector which appear 
most compatible with experiments at doping 1/8^ The 
density-wave part Hdw is given by 

^DW = E $ c4, C ^ ■ - + C '~ • - + h - c - 



- s kia k+QlM 



$l(ka) = -e*4 ^cos(fca; + ^) - cos fc^ 5t, 



a/2 l + >/2 



(2) 



s = (37r/4,7r), and the $ 2 s are 
y. <£ c implements a hop- 



where Q\ (tt/2,0), 
obtained from <I>J s by x 
ping modulation on the bonds of strength St, resulting 
in so-called valence-bond stripes j 13 i 14 Such bond-charge 
modulation, Fig. [TJi, is compatible with the STM data 
of Ref. and implies a strong d-wave component in the 
form factor F c of the CDW order parameter. The asso- 
ciated collinear spin order, Fig. [TJd, is implemented via a 
spin-dependent chemical potential Sfi a in <£ s ^ 

The pairing part Hp is dominated by a uniform d x i_ y i- 
wave pairing mean field Ao. In addition, a pairing modu- 
lation of amplitude SA is assumed along with the CDW, 
resulting in a pattern of bond pairing amplitudes qualita- 
tively similar to Fig. but with d-wave sign structure. 

k 

+ $l*ct c f ^ -fh.c., 
A(fc) = Ao (cos k x — cos k y ) , 
*£(fc) = -e^ i (cosfe + J) + cos ky) 5 A, (3) 

and & 2 is obtained from by x <->• 



III. INTER-LAYER JOSEPHSON COUPLING 

A. Inter-layer tunneling 

The Hamiltonian of the full system of two adjacent 
Cu02 layers is given by 



n = n 1 + n 2 + J2 \t±(k)cl ia c~ k2a + h.c. 



(4) 



ka 



with H l given in Eq. (pQ), and the inter-layer hopping 
matrix element^ 



t±(k) = ^( cos kx — cos k y ) 2 . 



(5) 



To calculate the Josephson coupling it is convenient to 
multiply global phase factors 6i to the superconducting 
mean fields Ao and SA in layer i. Then, the Josephson 
coupling measures the inter-plane phase stiffness: 

J j = \ [F(S0 = tt) - F(S6 = 0)] (6) 

where F = F 1 + F 2 + SF(S0), F % is the free energy of the 
isolated layer z, SF is the inter-layer tunneling contribu- 
tion to the free energy, and SO = #2 — 6\ . 

Assuming tj_ <C £, can be determined in second- 
order perturbation theory in t± : 



(7) 



a,/3=0 



where is the full Green's operator on layer z, T is the 
inter-layer tunneling operator from t±, f3 is the inverse 
temperature, and N the number of unit cells. The in- 
dices a,/? denote particle-hole space. Since the stripe 
directions are orthogonal, we have to consider the full 
(non-reduced) Brillouin zone (BZ). For a period-4 CDW 
(period-8 CDW+SDW) in each superconducting layer, 
the calculation of SF involves the diagonalization of a 
8x8 (16x16) Hamiltonian matrix to construct the 
Green's functions required in Eq. Q. 



B. Results 

The numerical calculations have been performed at 
zero temperature, with parameters t = 0.15 eV, Ao = 
0.024 eV, and fixed doping x = 1/8. For the homoge- 
neous case, this corresponds to \i = — 0.126 eV. 

Results for Jj for charge-only stripes as function of 
the hopping modulation St (SA = Sfi" 7 = 0) are shown 
in Fig. [2k • For large modulation amplitude, the Joseph- 
son coupling is seen to be strongly suppressed, e.g. by 
roughly a factor 10 for St = 0.07 eV. A simultaneous 
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FIG. 2: Inter-layer Josephson coupling Jj (a,b) and su- 
perconducting condensate amplitudes ipo (solid) and ip^ c 
(dashed) (c,d) as functions of the hopping modulation 
strength St, at fixed doping x — 1/8 and in the absence 
of magnetic order. The modulation in the pairing field is 
SA/St = (a,c) and 1.3 (b,d). The couplings and condensates 
are normalized w.r.t. the values of Jj and ipo at St — SA = 0. 



modulation in the condensate mean field by SA (here 
chosen to produce similar relative modulation strengths 
in the resulting bond kinetic energies and pairings) sup- 
presses the Josephson coupling even further, Fig. [2)3. 

Figs. [2]3 and d show the corresponding evolution of 
the homogeneous and modulated condensate amplitudes, 
ipo and ^q c , calculated from the solution of the mean- 
field Hamiltonian (pp). Here, we define ip as the sum 
of the magnitudes of the s-wave (on-site), s x 2 +y 2 -wave 
and d £C 2_ 2/ 2-wave condensates calculated from the real- 
space pairing amplitudes extracted from the solution of 
H l (note that the s and d x i_ y i representations of the 
point group mix in the presence of stripe order). A mod- 
ulated condensate ip^ c is always present for St ^ 0, but 
remains small if S A = 0. In contrast, for S A oc St as 
in Fig. [2]i, ip^ c increases and eventually dominates over 
ipo. A comparison between the evolution of Jj and ipo 
reveals that in the range of St where Jj drops dramati- 
cally, the uniform condensate ipo displays a much weaker 
depletion. This is also true for panels b and d where, at 
St « 0.08 eV, J j is reduced by a factor 200 while ipo still 
has half of its original value. From this we conclude that 
the primary source of the suppression of the Josephson 
coupling in our calculation is different from that of the 
PDW proposal by Berg et ali^ where the layer decou- 
pling is due to the absence of a homogeneous condensate. 

Analyzing our results further, we identify the 
momentum-space mismatch of the Fermi surfaces of the 
two layers, arising from the orthogonal stripe modula- 
tion, as the main source of the suppression of Jj. This 
mismatch is also accompanied by a mismatch of the 



nodal lines of the superconducting order parameter in 
the two layers, due to broken rotational symmetry in each 
layer. These effects can be nicely seen in the momentum- 
resolved contributions Jj(k) = SF^(S0 = tt) — SF^(S0 = 0) 
to the Josephson coupling, shown in Fig. [3j In the homo- 
geneous case, the largest contributions to Jj arise near 
the antinodal points of the (bare) Fermi surface. These 
contributions are drastically reduced (note the logarith- 
mic intensity scale) with increasing stripe modulation, as 
a result of the Fermi-surface distortions^^ accompany- 
ing the stripe order. For sizeable SA, the combination of 
Fermi-surface and order-parameter reconstruction even 
generates regions in momentum space with Jj(k) < 0, 
Fig. [3)3. (For S A = 0, this effect occurs only near the BZ 
diagonals due to the shift of nodal lines from stripe order, 
but this has little influence on Jj due to the specific mo- 
mentum dependence of the inter-layer tunneling.) The 
imposed pairing modulation SA is seen to contribute to 
the reduction of Jj, Fig. [4^. 

We now turn to the influence of magnetic SDW order 
as in Fig. [TJd. As shown in Fig. [4)3, SDW order alone 
(with CDW being parasitic only) leads to a moderate 
suppression of Jj. Similarly, SDW order in combination 
with a CDW suppresses Jj further compared to the non- 
magnetic case, mainly because of the additional Fermi- 
surface reconstruction arising from the SDW wavevector. 
Note that the relative spin orientation between the two 
layers does not enter the result. 

Finally, we link our findings to the experimental situ- 
ation. Unfortunately, the magnitude of the charge mod- 
ulation in cuprate stripes is not well known: From res- 
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FIG. 3: Positive (top) and negative (bottom) momentum- 
resolved contributions to Jj. Each panel shows \Jj(k)/t\\ 
(see text) as function of k on a logarithmic intensity scale. 
The modulation strength St is zero in a) and increases from 
b) to c). The other parameters are as in Fig. [2j),d. From 
a) to c) the positive contributions near the antinodal points 
are reduced, moreover negative contributions appear. Note 
that the momentum dependence of the inter-layer tunneling 
(0 suppresses the contributions near the diagonals. 
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FIG. 4: Normalized inter-layer Josephson coupling Jj at dop- 
ing 1/8 for a) paramagnetic stripes with hopping modulation 
St and accompanying pairing modulations SA/St = (solid), 
0.5 (dashed), 1.3 (dotted), and b) antiferromagnetic stripes 
with magnetic modulation set by for St = 5 A = (solid) 
and for moderate fixed kinetic energy modulation of ~ ±25% 
(dashed) achieved via adjusting St and keeping SA = 1.3St. 

onant soft x-ray scattering^ on La 15 / 8 Ba 1 / 8 Cu04 the 
modulation on the oxygens was concluded to be of order 
factor 4, but the quantitative analysis is model depen- 
dent. From the STM data 7 on Bi2Sr2CaCu208+5 and 
Ca2-^Na ;E Cu02Cl2 one may infer a typical modulation 
amplitude in the charge sector of ±20... 30%. In our 
calculation, a modulation of ±25% in the bond kinetic 
energies is obtained from St = 0.023 eV, which gives a 
20% reduction of Jj (Fig.[2)o), while St = 0.07 eV with a 
factor 90 reduction of Jj corresponds to a kinetic-energy 
modulation of about a factor 10 (here, the neglect of 
Mott physics makes this number less meaningful, as the 
quasiparticle theory does not account for half filling being 
special). 

The magnitude of the magnetic order in striped 214 
cuprates is known reasonably well, at x = 1/8 the max- 
imum moment size is 50 . . . 60% of that of the undoped 
compound^^ In our mean-field calculation, Sfi a « 
±0.07 eV corresponds to a maximum moment of 0.36/i# 
and gives an additional factor 2 suppression of Jj. 



As stripes are particularly stable near doping 1/8, it 
is conceivable that Lai5/ 8 Ba!/ 8 Cu04 has rather strong 
modulation in the charge sector, in which case Jj would 
be suppressed drastically by stripe order. This in turn 
would be consistent with the fact that the 2d fluctuating 
pairing regime is restricted to the vicinity of x = 1/8. Of 
course, corrections beyond the mean-field picture arising 
from phase fluctuations, quenched disorder, short-range 
magnetism, and other pseudogap physics is of impor- 
tance to fully understand the hierarchy of energy scales 
in La 15 / 8 Ba 1 / 8 Cu04. 



IV. CONCLUSIONS 

Using a mean- field quasiparticle framework, we have 
calculated the inter-layer Josephson coupling in super- 
conducting stripe states. We have assumed that the in- 
plane stripe orientations alternate from layer to layer, 
as is the case in 214 cuprates with LTT lattice structure. 
For realistic stripe modulation strengths, we find that the 
inter-layer coupling can be easily reduced by an order of 
magnitude. The primary cause of this reduction is the 
momentum-space mismatch between the reconstructed 
Fermi surfaces of adjacent layers, while the depletion 
of the zero-momentum superconducting condensate (in 
favor of a modulated one) is secondary. Whether this 
effect is sufficient to explain the unusual properties of 
La 15 / 8 Ba 1 / 8 Cu04 is not yet clear. 
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